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Представлены аналитические выражения и результаты численных расчетов углового и частотного спектров радиосигнала, рассеянно- 
го на внешней поверхности турбулентного плазменного образования, создаваемого реактивным двигателем космического аппарата. 


Введение 

В работах [1,2] предложена модель турбулентно- 
го плазменного образования, создаваемого реактив- 
ным двигателем космического аппарата (КА) в ио- 
носфере. Полагается, что турбулентные плазменные 
неоднородности сосредоточены вблизи поверхности 
тел вращения: конуса, параболоида вращения и по- 
верхности, образованной вращением кривой четвер- 
того порядка. Рассеяние плоской радиоволны на 
внутренней поверхности полого турбулентного плаз- 
менного образования в случае радиозондирования 
вслед факелу ракеты, исследовано в работах [2-9]. 

Представляет интерес исследование рассеяния 
радиоволн КВ и УКВ диапазонов на внешней пове- 
рхности турбулентного плазменного образования, 
создаваемого реактивным двигателем КА в ионос- 
фере, результаты которого излагаются ниже. Рас- 
сматривается рассеяние плоской радиоволны на 
турбулентных плазменных неоднородностях, сос- 
редоточенных вблизи поверхности тел вращения: 
конуса и параболоида [2, 9-14]. 

Предполагается, что размеры тел вращения 
намного превосходят длину радиоволны и харак- 
терные размеры турбулентных неоднородностей, а 
час тота р адиоволны удовлетворяет соотношению 
/ <гЩр%/ р где /V]. и А 0; - возмущенная и невозму- 
щенная концентрации заряженных частиц ионос- 
ферной плазмы, / р - плазменная частота ионосферы 
на высоте движения КА. При этом, как отмечалось 
в работах [1-3], коэффициент отражения от турбу- 
лентных плазменных неоднородностей близок к 1. 

Полагается, что турбулентные неоднородности 
распределены равномерно по поверхности тел вра- 
щения. Они имеют как направленную скорость, 
ориентированную вдоль образующей тел вращения, 
так и случайную изотропную скорость, распреде- 
ленную по нормальному закону. Исследуются три 
типа диаграмм переизлучения турбулентных неод- 
нородностей: изотропная, ламбертовская и квази- 
зеркальная. Решение задачи проводится на основе 
применения теории переноса излучения [15]. 

1. Частотный спектр сигнала, рассеянного 

на внешней поверхности турбулентного 

плазменного тела вращения 

Для решения задачи применяется сферическая 
система координат (рис. 1), центр которой распо- 
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ложен в вершине тела вращения, а ось х направле- 
на вдоль его оси, Ѳ - полярный угол, (р - азимуталь- 
ный угол, отсчитываемый от плоскости, имеющей 
направление распространения падающей волны е, 
и ось х . Индекс и соответствует координатам пове- 
рхности тела вращения, і - параметрам падающей 
волны, а — рассеянной. 



Рис. 1. Система координат, используемая при расчетах рас- 
сеяния радиоволн на турбулентном теле вращения 

Диаграмма переизлучения турбулентностей, не 
зависящая от направления распространения пада- 
ющей волны е, , может быть представлена в виде [2] 

р Ю = А ] (ё $ яу, (і) 

где е [ - единичный вектор в направлении распрост- 
ранения рассеянной волны, л - нормаль к поверх- 
ности, А і - коэффициент, который вычисляется из 
условия нормировки \Р(е)сК1=\, сІО. - элемент те- 
лесного угла, в который происходит рассеяние. Из 
условия нормировки следует А/=(]+\)/2л. Случай 
7=0 соответствует изотропной диаграмме переизлу- 
чения турбулентностей, а /= 1 - ламбертовской. 

Представляет интерес рассмотрение диаграммы 
переизлучения турбулентностей квазизеркального 
типа, учитывающей направление распространения 
падающей волны ер. 

Р(ё,) = А( Ае 1 2 т -Ае 2 ), (2) 

где Ае = - ё Е0 , ё х0 = е. - 2 Я (Яё, ). 

В этом случае диаграмма переизлучения асим- 
метрична, основная часть энергии рассеянного 
сигнала сосредоточена в окрестности вектора е 
который соответствует направлению зеркального 
отражения волны от поверхности. Выражения для 
Ае 2 , Ае 2 и А в (2) имеют следующий вид: 

Ае 2 = 2(1-еД +2(пё і )(пё і )), (3) 
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\е 2 т = 2(\ + (\-(пё і ) 2 ) ш ), (4) 

А = і/[4тг(^/і — (ле,-) 2 -щ)\, (5) 

Выражение для спектральной плотности энер- 
гии, рассеянной на турбулентном плазменном теле 
вращения, следуя [14, 15], можно записать в виде 

сП / 1 0 = -(ё і п )Р(ё, )(1ё>(ПИ 2л, ( 6 ) 

где сП, - количество энергии, рассеянной в задан- 
ном направлении, / 0 - энергия падающей волны, 
сі В - элемент рассеивающей поверхности. 

При рассеянии электромагнитной волны на пе- 
ремещающихся турбулентностях происходит допле- 
ровский сдвиг частоты А/ Удобно ввести величину 
безразмерного доплеровского сдвига частоты рассе- 
янного сигнала, которая определяется выражением 

/. = (Мс)ИШ = 

= (ё,-ёма+з ѵ/ѵ„)-цѵ 0 /ѵ и ), ( 7 ) 

где с - скорость распространения электромагнит- 
ной волны в плазме, То - частота падающей волны, 
Е 0 - скорость перемещения КА, Ѵ„ и 5Ѵ- направ- 
ленная скорость и скорость случайного перемеще- 
ния турбулентностей вдоль образующей тела вра- 
щения, е и и е, единичные векторы, направленные 
вдоль образующей и оси потока соответственно [2]. 

В сферической системе координат выражение 
(7) записывается в виде 

/„ = (8ІП Ѳ, С08 ср, - 8ІП Ѳ і ) 8ІП Ѳ С08 (р„ + 

+ 8ІП Ѳ, 8ІП (р 5 8ІП Ѳ и 8ІП (р и + 


/» = 2((У 0 / У и ) С 08 Ѳ і -(1 + ЗѴ/Ѵ и )х 
х(зіп Ѳ, С 08 (р и + 2а,р С 08 Ѳ і ) / 0, 


где <2 = Л /і + 4о 1 р 2 . 

Определяя из (10) зависимость ср и (/,) и подстав- 
ляя ее в (9), получаем выражение для частотного 
спектра рассеянного назад сигнала: 

я/:і = - ад + в) 2 у т , ( і о 

ѵОЬ КУ и 


где О 
В 


1 


2(1 + <5 У / У и ) 8Іп Ѳ і зіп Ѳ и ’ 


1 


У и 8ІпѲ і 1 ё ѲД1 + «5 Ѵ/Ѵ и ) і ё Ѳл ё Ѳ ’ 
ё,ёі = 8ІП Ѳ ; С08 Ѳ и (I)/» + В) - С08 6) 8ІП Ѳ и . 


Выражение (11) определяет текущий частотный 
спектр рассеянного радиосигнала. Однако при экс- 
периментальных измерениях энергия принимаемого 
сигнала для заданной частоты / определяется в час- 
тотном интервале А/>0, поскольку согласно соотно- 
шению неопределенностей А/Аіа 1 , а время анализа 
А/ составляет конечную величину (Д/«10...100 с). 

Для получения частотного спектра необходимо 
проинтегрировать по / выражение (11) в заданном 
частотном интервале от Ц до Ц+А/, и поделить на 
величину этого интервала 


Ш) = о 


А/, С08 Ѳд 


Л.+У. 

1 


Л. 


(е,п)Р(ё,)с1/, / 12) 

(1 -іур + В )У 2 ' 


+(С 08 Ѳ, — С 08 Ѳ,-)(со8 Ѳ і - Ѵд/ Ѵ и ). (8) 

Решая ур. (8) относительно азимутального угла 
ср и и подставляя полученную зависимость ср и (/,) в 
ур. (6), можно получить выражение для частотного 
спектра радиоволн, рассеянных на плазменном об- 
разовании 


Для диаграмм переизлучения рассеивателей по 
закону Ламберта и изотропного типа интеграл в 
правой части ( 12) вычисляется аналитически, и вы- 
ражение (12) принимает вид 

^ ш = в Ѳц тл, +ю- щ/ к . и, (із) 


Ш) = 


сП, Ал 


{ §ѳ п 

С08 Д 


Щ) 


ёсрЛЛ) 


(9) 


где $„(/,) - нормированная величина спектральной 
плотности энергии, рассеянной в единичный эле- 
мент телесного угла, і т - размер тела вращения 
вдоль оси і. 

В случае рассеяния назад (е=-е,) для коничес- 
кого потока из (8) следует 

/,=2((У 0 /^)со8Ѳ,-(1 + 5У/У„)х 
х(с 08 Ѳ, С 08 Ѳ и + 8ІП Ѳ, 8ІП Ѳ„ С 08 (р„ )), (10) 

где Ѳ и - угол полураскрыва конуса, ср и - азимуталь- 
ный угол, определяющий положение турбулент- 
ности на поверхности потока, Ѳ, - угол между нап- 
равлением распространения падающей волны ё, и 
осью і потока. 

Если поток имеет форму параболоида враще- 
ния і=а { р г , то 


где при /=0 

В>Ш = 5ІП Ѳ, С08 Ѳ и (ВО, -О) - 
-О г (В 8ІП Ѳ, С08 Ѳ і - С08 Ѳ і 8ІП Ѳ и ), '^) 

Д = агс8Іп(-Е)/, - В), (15) 

й = л/і ~(/*0 + В) 2 . (16) 

При 7=1 величина (Ѵв (13) имеет вид 

ЩУ) = (/,О~Щ0/2 + 

+2 (В 2 + 1)Д зіп 2 Ѳ, С 08 2 Ѳ„ + (ВЦ. - О) X 
х(258Іп 2 Ѳ С 08 2 Ѳ -8Іп2Ѳ. 8Іп2Ѳ /2) + 

к і и і 

+Д, (В 2 8ІП 2 Ѳ і С08 2 Ѳ и + 

+ С 08 2 Ѳ і зіп 2 Ѳ и - В 8Іп261 зіп 2Д /2). 

В случае диаграммы переизлучения квазизер- 
кального типа Р(е ,) определяется выражением (2) и 
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интеграл в правой части (12) аналитически не вы- 
числяется, вследствие чего проводится численное 
интегрирование методом Монте-Карло. 

На рис. 2 приведены расчеты зависимости час- 
тотного спектра сигнала, рассеянного на коничес- 
ком турбулентном потоке, которые проводились для 
значений параметров Ѵ й /Ѵ=2, 0 О =26,6°, 5Ѵ/Ѵ= 0. 

Сплошные кривые на рис. 2 рассчитывались по 
формуле (11) для различных значений угла зондиро- 
вания потока Ѳ,. Гистограммы соответствуют частот- 
ному спектру, вычисленному по формулам (12-17). 


Частотный спектр рассеянного сигнала харак- 
теризуется монотонным возрастанием энергии с 
увеличением безразмерного доплеровского сдвига 
частоты до значения 


/, = 2 


г^-созѲ, -(1 + 8Ѵ /Ѵ и )со&(Ѳ. +Ѳ и ) 

V Ѵ и 


(18) 


Выражение (18) следует из (8) при <р= л, а при 
(р=л/2, 3 лг/2 наблюдается минимальная энергия 
рассеянного назад сигнала. Соответствующее вы- 
ражение для /. имеет вид 


(у 

/„ =2 — ^-созѲ ~(\+8Ѵ ! С )со8 0 соз 0 

V 'и 


(19) 



Рис. 2. Частотный спектр радиосигнала, рассеянного на 
внешней поверхности турбулентного плазменного 
образования 


На рис. 3 приведены результаты расчета верхних 
и нижних частот рассеянного сигнала, формулы 
(18, 19). По горизонтальной оси отсчитан угол 0, 
между осью і конического потока и направлением 
падающей волны, а по вертикальной оси - значе- 
ние безразмерного доплеровского сдвига частоты /,'. 

Сплошные кривые соответствуют отношению 
скорости перемещения источника к скорости рас- 
сеивателей потока Ѵ й /Ѵ=2, штриховые - Ѵ 0 Ѵ=2,5 и 


штрихпунктирные - Ѵ 0 /Ѵ=3. Величина доплеровс- 
кого сдвига частоты уменьшается с ростом угла 
зондирования потока 0, и увеличивается с ростом 
отношения скоростей Ѵ п /Ѵ и . Расчеты проведены 
для значений угла полураскрыва конуса 0„=26,6°. 



Рис. 3. Зависимость верхних и нижних частот рассеянного 
радиосигнала от полярного угла 0 ( и отношения ско- 
ростей Ѵ,/Ѵ„ 

Из выражений (18, 19) следует соотношение для 
ширины полосы частот сигнала, рассеянного на 
коническом потоке: 

/ ве =2$ + 8Ѵ/Ѵ а )зтѲ і зтѲ я . (20) 

Из (20) можно получить соотношение для ши- 
рины полосы частот сигнала, рассеянного на пото- 
ке, имеющем форму параболоида вращения і=а { р 2 

/в Р =2(\ + 8Ѵ/ Ѵ и ) 8Іп Ѳ, / л /і+4а 1 2 р 2 - 

Отношение / 8с І/ 8р = 8Іп0„ ^1 + 4а;р 2 . 

Так, при малых значениях координаты располо- 
жения рассеивателя на поверхности параболоида 
р<1/(2а 1 іё,Ѳ и ), величина /&.</%, а с увеличением р 
выполняются обратные соотношения, и ширина 
полосы частот сигнала, рассеянного на параболои- 
де, уменьшается. 


2. Сечение рассеяния радиоволн на турбулентном 
плазменном конусе 


Выражение для дифференциального эффектив- 
ного сечения рассеяния (имеется в виду рассеяние 
в элемент телесного угла сЮ. в заданном направле- 
нии е,) электромагнитной волны на турбулентном 
плазменном образовании, согласно (6), имеет вид 

^ = -\Ж(ёр ) Р(ё і ), ( 21 ) 


Для конической поверхности (21) преобразует- 
ся в выражение 


сіа _ 
сЮ. 


2 2 1§0 

т и 

2 со8 Ѳ и 


|г/<р„(ей)Р(еД, 


а 


( 22 ) 


где і т размер конуса вдоль оси і. Пределы интегри- 
рования по (р и в (22) определяются из совместного 
решения системы неравенств 


ер < 0, 

(23) 

еТі > 0, 

(24) 
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первое из которых является условием облучения 
данного элемента поверхности падающей волной, 
а второе - условием облучения данным элементом 
рассеивающей поверхности приемника. 

В результате подстановки ( 1 ) в (22) для изотроп- 
ной по полусфере диаграммы переизлучения тур- 
булентностей (/'=0) получаем выражение для сече- 
ния рассеяния: 

сіа _ 2 т 1 8@„ х 
сІО. 4 л 


х[созѲт ё Ѳ X <Р„ |!і -8ІпѲ„ X 8ІП % |а‘]- (25) 

к= 1,2 і=1,2 


В случае у— 1 в формуле (1) рассеяние происхо- 
дит в соответствии с ламбертовской диаграммой 
переизлучения турбулентностей и выражение для 
сечения рассеяния принимает вид 


сіа 

ап 


2 -А [в ^ 5 - т2(Ри 


к= 1,2 


+С X С 08 2(р„ +Б X 8Іп<р„ |*‘ + 


+Е ^ С 08 <р„ ь а к к +Е^(р„ Ъі 


(26) 


где 5 = - С08 Ѳ и 8 Іп Ѳ і 8 Іп Ѳ,. С08 ( р х / 2, 

С = С 08 8 ІП О, 8 ІП Ѳ Л 8 ІП ТД / 2, 

I) = 2 8ІП (С08 Ѳ і 8ІП Ѳ К С08 (р х +8ІП 6* С08 6( ), 
Іі=-28ІпѲ С08 Ѳ 8ІП Ѳ 8ІПЮ , 
і 7 = -28ІП0 1 ё 0 С08 Ѳ 8Іп Ѳ +2 В. 

При подстановке (2-5) в (22) получаем выраже- 
ние для сечения рассеяния в случае закона переиз- 
лучения турбулентностей квазизеркальнош типа: 


сіа _ 
сЮ. 




87Г С08 Ѳ, 


г 

| гіфи (Щп) 


20+ф-(ё і п) 2 )-Ае 1 

х/!-(е,«) 2 -Ц 


(27) 


нат, ср, - полярный угол, Ѳ . ; - длина радиус-вектора. 
В левом столбце на рис. 4-6 представлен вид спере- 
ди п/2<Ѳ<п, а в правом - вид сзади 0<Ѳ<я/2. 

Рис. 4 получен для изотропной по полусфере ди- 
аграммы переизлучения шероховатой поверхности 
7=0. С ростом О, проекция поверхности конуса, об- 
лученной падающей волной, на плоскость попереч- 
ную ё і вначале увеличивается, а затем, начиная с 
Ѳ,>Ѳ" ( Ѳ * зависит от угла полураскрыва конуса 0„), 
уменьшается, что приводит к аналогичной зависи- 
мости величины сечения рассеяния от угла Ѳ 



Рис. 4. Индикатриса рассеяния радиоволн на турбулентном 
конусе в случае изотропной диаграммы переизлуче- 
ния турбулентностей 



Рис. 5. Индикатриса рассеяния радиоволн на турбулентном 
конусе в случае диаграммы переизлучения турбуле- 
нтностей по закону Ламберта 


Вычисление интеграла в правой части (27) про- 
ведено методом Монте-Карло. При вычислении 
учитывается, что вклад в сечение рассеяния дает 
лишь область рассеивающей поверхности, облу- 
ченная падающей волной и облучающая прием- 
ник, определяемая из совместного решения систе- 
мы неравенств (23, 24). Относительная погреш- 
ность вычислений составляла величину <0,05. 

Результаты численных расчетов индикатрисы 
рассеяния электромагнитных волн на шероховатой 
поверхности конуса для рассмотренных типов ди- 
аграмм переизлучения шероховатой поверхности 
(25-27) представлены на рис. 4-6 в виде изолиний 
сісг/сЮ. 
величины — - — . 

Расчеты проведены для значений угла полураск- 
рыва конуса Ѳ„=20°. Вверху в центре на рис. 4-6 по- 
казаны соответствующие диаграммы переизлучения 
элементов шероховатой поверхности. Результаты 
расчетов представлены в полярной системе коорди- 


Область поверхности, облучающая приемник, 
и, соответственно, величина сечения рассеяния, 
увеличиваются с ростом угла рассеяния Ѳ \ и дости- 
гают максимума при Ѳ>л-Ѳ„, после чего остаются 
постоянными до Ѳ=л. Постоянство сечения рассе- 
яния при л-Ѳ<Ѳ< л объясняется тем, что в указан- 
ном случае вся поверхность, облученная падающей 
волной (е,п <0), облучает приемник (е,п>0). 

В случае ламбертовской диаграммы переизлуче- 
ния турбулентных неоднородностей (рис. 5) направ- 
ление максимума индикатрисы определяется нап- 
равлением нормали к поверхности, в окрестности 
которой сосредоточена энергия излучения. При 
этом максимум индикатрисы реализуется при следу- 
ющих значениях углов: Ѳ-=л/2-Ѳ и , Ѳ=л-Ѳ-=л/ 2+Ѳ„. 

Диаграмма переизлучения турбулентностей ква- 
зизеркального типа (рис. 6) обладает максимумом, 
ориентированным в направлении зеркального отра- 
жения волны от поверхности. При этом максимум 
индикатрисы располагается под углом 0/=0 ( + 2Ѳ„. 
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Рис. 6. Индикатриса рассеяния радиоволн на турбулентном 
конусе в случае квазизеркальной диаграммы переиз- 
лучения турбулентностей 


Для всех рассмотренных диаграмм переизлуче- 
ния с увеличением угла тр, от 0 до л - величина сече- 
ния рассеяния растет и достигает максимума при 
(р=л, так как лишь в этом случае все элементы рас- 
сеивающей поверхности, расположение которых 
удовлетворяет условиям (23, 24), вносят свой вклад 
в рассеянный сигнал. 

Расчеты, проведенные для других значений угла 
полураскрыва конуса Ѳ и , показывают, что величина 
сечения рассеяния возрастает с ростом угла Ѳ ю что 
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объясняется увеличением освещенной волной об- 
ласти поверхности, участвующей в формировании 
рассеянного сигнала. 

Заключение 

Исследовано рассеяние радиоволн на турбулент- 
ных неоднородностях слабоионизованной плазмы, 
сосредоточенных вблизи поверхности тел враще- 
ния: конуса и параболоида. Рассмотрены три типа 
диаграмм переизлучения турбулентных неоднород- 
ностей: изотропная, ламбертовская и квазизеркаль- 
ная. Получены аналитические выражения и прове- 
дены численные расчеты частотного спектра и сече- 
ния рассеяния радиоволн на расширяющемся тур- 
булентном потоке. При возвратном зондировании 
частотный спектр сигнала, рассеянного на коничес- 
ком потоке, характеризуется монотонным возраста- 
нием энергии с увеличением доплеровского сдвига 
частоты. Величина доплеровского сдвига частоты 
уменьшается с ростом угла зондирования и с умень- 
шением отношения скорости движения КА к нап- 
равленной скорости перемещения турбулентных 
неоднородностей. Ширина полосы частот рассеян- 
ного сигнала возрастает с увеличением угла зонди- 
рования потока и относительной дисперсии скорос- 
ти перемещения турбулентных неоднородностей. 
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